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Abstract 

Lattice models with long-range interactions of power-law type are suggested as a new 
type of microscopic model for fractional non-local elasticity. Using the transform oper¬ 
ation, we map the lattice equations into continuum equation with Riesz derivatives of 
non-integer orders. The continuum equations that are obtained from the lattice model 
describe fractional generalization of non-local elasticity models. Particular solutions and 
correspondent asymptotic of the fractional differential equations for displacement fields 
are suggested for the static case. 

PACS: 45.10.Hj; 61.50.Ah; 62.20.Dc 


1 Introduction 

Lattice with long-range interaction is a subject of investigations in different areas of mechanics 
and physics (see for example [D El El la [5] and PEI El E]). The long-range interactions have 
been studied in discrete systems as well as in their continuous analogues. As it was shown in 
mm (see also mu), the continuum equations with derivatives of non-integer orders can 
be directly connected to lattice models with long-range interactions of power law type. 

The theory of derivatives and integrals of non-integer orders mm allow us to investigate 
the behavior of materials and media that are characterized by non-locality of power-law type. 
Fractional calculus has a wide application in mechanics and physics (for example see [T51 ITBl 
EZl EHl n El OSl ED]). The fractional calculus allows us to formulate a fractional generalization 
of non-local elasticity models in two forms: the fractional gradient elasticity models (weak 
power-law non-locality) and the fractional integral non-local models (strong power-law non- 
locality). Fractional models of non-local elasticity and some microscopic models are considered 
in different articles (see for example [211 E21 E31 EH ESI ES], [23 ESI EH ED] and [211 E21 E3]) 
Elastic waves in nonlocal continua modelled by a fractional calculus approach are considered in 
[?, EH 123 ES]- In [23 E2] a general approach to describe lattice model with power-law spatial 
dispersion for fractional elasticity has been proposed. This approach can be used for different 
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type of interaction of lattice particles. Therefore explicit forms of the long-range interactions 
are not considered in [3ll|32]. In |33] a model of lattice with long-range interaction of Griinwald- 
Letnikov-Riesz type has been snggested to describe fractional gradient and integral elasticity 
of continnnm. In this paper we focns on the lattice models with long-range interaction of 
power-law type as new type of microscopic models for fractional generalization of elasticity 
theory. We snggest lattice models with power-law long-range interaction as microscopic model 
of fractional non-local continnnm. The eqnations for displacement held are directly derived from 
the snggested lattice models by the methods of mm- The snggested generalization of the 
elasticity eqnations contains the fractional Laplacian in the Riesz’s form [H]. We demonstrate 
a connection between the dynamics of lattice system of particles with long-range interactions 
and the fractional continnnm eqnations by using the transform operation suggested in mm- 
We show how the continuous limit for the lattice with long-range interactions of power-law 
type gives the continuum equation of the fractional elasticity. We get particular solutions of 
the fractional differential elasticity equations for some special cases. 


2 Equations of lattice model 


As a microscopic model, we use unbounded homogeneous lattices, such that all particles are 
displaced from its equilibrium positions in one direction, and the displacement of particle is 
described by a scalar held. We consider one-dimensional lattice system of interacting particles. 
The equations of motion for particles are 




+ CO +00 

92 ^ K 2 {n,m) [un- Ur^ + Qa ^ Ka{n, m) {un - Ur^ + F{u), (1) 

m=—00 m=—00 

m^n m^n 


where Un(t) are displacements from the equilibrium, g 2 and Qa are the coupling constants of 
particle interactions, and the terms F{n) characterize an interaction of the particles with the 
external on-site force. For simplicity, we assume that all particles have the same mass M. The 
function K 2 {n,m) describes the nearest-neighbor interaction with coupling constant 92 = K, 
which is the spring stihness. The function Ko,{n,m) describes the long-range interaction with 
a coupling constant 9 o,- For a simple case each particle can be considered an inversion center 
and 

Ka{n, m) = Ka{\n — m\). 

Equations of motion ([T]) have the invariance with respect to its displacement of lattice as 
a whole in case of absence of external forces. It should be noted that the non-invariant terms 
lead to the divergences in the continuous limit |5]. 

Using the approach suggested in [UllinilS], we can consider a set operations that transforms 
the lattice equations for Unit) into continuum equation for displacement held M(a;, t). We assume 
that Unit) are Fourier coefficients of the held uik,t) on [—/co/2, A:o/2] that is described by the 
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equations 


1 |•+ko/2 

Unit) = Tr u{k,t) = J'^^{u{k,t)}, 

ko J-ko/2 

-\-oo 

u{k,t) = ^ Unit) = J^A{Un(t)}, 


( 2 ) 

(3) 


where Xn = nd and d = 2n/ko is distance between equilibrium positions of the lattice particles. 
Equations (jSj) and ([2]) are the basis for the Fourier series transform and the inverse Fourier 
series transform J-^^. 

The Fourier transform can be derived from ([2]) and ([2]) in the limit as d —?■ 0 {ko —>■ cxo). In 
this limit the sum is transformed into an integral, and equations ([2]) and ([3]) become 


u{k, t) 




dx e ^^^u{x,t) = T{u{x,t)}, 


(4) 


-I ^+oo 

u{x,t) = — / dk e^^^u{k,t) = T~^{u{k,t)}. (5) 

27 r J_oo 

Here we use the lattice function 

/ \ 27 r - 

Un{t) = -^u{Xn,t) 

Kq 

with continuous function u{x,t), where Xn = nd = {27in)/ko —)■ x. We assume that u{k,t) = 
Cu{k,t), where C denotes the passage to the limit d —)■ 0 (fco ^ oo), i.e. the function u{k,t) 
can be derived from u{k,t) in the limit d —)■ 0. Note that u{k,t) is a Fourier transform of 
the held u{x,t). The function u{k,t) is a Fourier series transform of Un(t), where we can use 
Unit) = {27i/ko)u{nd,t). 

We can state that a lattice model transforms into continuum model by the combination 
C J-A of the following operation mm-- 
The Fourier series transform: 


J^A ■ Unit) -> J^A{Un(t)} = u{k,t). 

The passage to the limit d ^ 0: 


C = lim : u{k,t) C{u{k,t)} = u{k,t). 

d-^O 

The inverse Fourier transform: 

u{k,t) ^ J^~^{u{k,t)} = u{x,t). 


( 6 ) 

(7) 

( 8 ) 


These operations allow us to get continuum equations from the lattice equations [lllDroiE]. 
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3 Lattice with nearest-neighbor interaction 

Let us consider the lattice with nearest-neighbor interaction that is described by ([1]), where 
Ka{n — m) = 0, and 


H-oo 

K2{n,m) Um{t) = Un+l{t) -2Un{t) +Un-l{t), (9) 

m=—oo 

m^n 

where the term K 2 {n,m) describes the nearest-neighbor interaction. Let us derive the usual 
elastic equation from the lattice model with the nearest-neighbor interaction with coupling 
constant g 2 = K, which is the spring stiffness. The following statement [luiiniis] gives for this 
lattice model with the nearest-neighbor interaction the corresponding continuum equation in 
the limit d —?• 0. 

Proposition 1. In the continuous limit (d ^ Q) the lattice equations of motion 

= K ■ (^nn+i(t) - 2nn(f) nn-i(^)) + F{n) (10) 

are transformed by the combination IF~^C of the operations (013) into the continuum equa¬ 
tion: 

^ ^g^2 ^U{X, t) + ^f{x), ( 11 ) 

where 

M E Kd"^ Kd 

and C‘1 is a finite parameter, A is the cross-section area of the medium, E is the Youngs mod¬ 
ulus, and f{x) = E[x)/{Ad) is the force density. 

A detailed proof of Proposition 1 is given in Appendix 1. 



As a result, we prove that lattice equations fITU]) in the limit d —)■ 0 give the continuum 
equation with the Laplacian (see also m)- Note that this result can be derived by methods 
described in Section 8 of |12], where the relation 


exp i 



u{x, t) 


u{x -\- d, t) 


and the representation of flTOl) by pseudo-differential equation are used. 
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4 Lattice model with long-range interaction 


Let us derive a continuum equation for the lattice with long-range interaction that is described 
by O, where K^iji — m) satishes the conditions 

OO 

Ka{n — m) = Ka{\n — m\), E \Ka{n)\'^ < OO. (13) 

n=l 

To have fractional gradient elasticity models, we assume that the function 


H-OO 


Ka{k) = e = 2 iLo,(u) cos(fcn), 


(14) 


n=—oo 

n^O 


n=l 


satishes the condition 




k^O |/c|“ 

where Aa has a hnite value. Condition flTHll means that 


(15) 


- /q(0) = A^\k\‘‘ + /?„(«:), 


(16) 


for A; —)• 0, where 

lim /?„(A))/|fc|“ = 0. 


(17) 


The interaction terms iLodn—m|), which give the continuum equations of gradient elasticity 
models, can be dehned as 


KJn) = 


-1)' 


r(Q;/2 -|- 1 -|- u)r(Q;/2 -|- 1 — n) 
Using the series (Ref. |13], Sec.5.4.8.12) 


(18) 


E 


-1)' 


■)2y-l 


^ r(z/ -|- 1 -|- n)V{v -|- 1 — u) 
where v > —1/2 and 0 < A; < 27r, we get 


cos(uA)) = 


r(2r7 -|- 1) 


sm 


2v 


2 ^““^ fk 

Kaik) — KJ^) = — -- sin“ ( — I = 

^ ^ ^ ’ r(a + l) \ 2 ) 2 r(a + l) 


2r2(z/ + 1)’ 


A:“ + 0(A:“+2). 


(19) 


( 20 ) 


Here we use z/ = a/2 and sin(A:/2) = k/2 + 0{k^). The limit A: —)■ 0 gives 


K^{k)-k^{0)_ 1 

fc™ |A)|“ 2r(a + l)’ 


( 21 ) 
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and we have Aa = l/{2T{a + 1)). To consider a fractional generalization of the elastic theory, 
the variables x and d = Ax are dimensionless. Note that the interaction flTSl) for integer valnes 
of a is discussed in [H]. 


Proposition 2. The lattice equations 
d'^u. 


+ 0O 


M- 


dt"^ 


= g 2 i^ 2 (n, m) [u^ - + 9a K^{n - m) [u^ - + F{n), (22) 


m=—oo 
m^n 


m=—oo 

m^n 


where g 2 and ga are coupling constants, K 2 {n, m) is defined by m and Ka{\n — m\) is defined by 
m. are transformed by the combination iF Fa of the operations into the continuum 
equation: 


^ = *^2 Au{x, t) - Ca {-A)°Fu{x, t) + -f{x), 

ott p 

where {—A)°‘F fractional Laplacian in the Riesz’s form and 

92 d^ „ gad"" 


C2 = 


AM ’ 


C'„ = 


2r(a + 1)M 


(23) 


(24) 


are finite parameters. 

A detailed proof of Proposition 2 is given in Appendix 2. 

In the Proposition 2, we use the Riesz fractional derivative (—A)“/^. It can be dehned as 
non-integer power of the Laplace operator in terms of the Fourier transform F by 

((-t^rig)(x) = F~d\K(rm) 


(25) 

This fractional Laplacian can be also dehned in the form of the hypersingular integral [I31[I1] 
by 


{{-AT^ffix) = 


:{ATf){z)dz, 


dn{m,a) jRn \zfi+^ 

where m > a > 0, and (A™/)(z) is the hnite difference of order m of a function f{x) with a 
vector step z G M"" and centered at the point x G M”: 




m\ 


k=0 


k\{m — k)\ 


f{x - kz), 


where the constant dn{m,a) is dehned by 

dn{rn, a) = 


F+^FAm{a) 


2"P(1 -|- a/2)T{n/2 + oi/2) sin(7rQ;/2) ’ 
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and 


i=o 


m\ 






This hypersingular integral does not depend on the choice of m > a > 0 . 

We can note another possibility to set the interactions described by K 2 {n,m). The term 
K 2 {n, m) that describes the nearest-neighbor interaction can be represented in forms that differ 
from Q. In general K 2 {n,m) describes a special form of the long-range interaction. Let us 
give some example of these forms. 

Example 1. Instead of the nearest-neighbor interaction function ([9]) we can use the long- 
range interaction with 


we obtain 


Then we have 



K2{n,i 

. 5.4.2.12) 


y - - — cos[nk) 

1 ^ 

n=l 

+00 

K2{k)^2j2 

n=l 

(-1)" 

2 ' 


Mk) 


\n — m\ 


cos{nk) = - { — 3 ” ^ ’ 1^1 ^ 


1 TT^ 

cos(fcn) = -/c^ ——, \k\ ^ 71. 

2 6 


Example 2. If we consider the long-range interaction in the form 

(_l)rx+l 


Kon) = 


m? — 


then equation fflTl) gives 

For —)■ 0, we obtain 


k{k) = —— -cos{ak) - 
asm vra 


k2{k)-k2{o) = 


Tia 


2sinf7ra' 


k + o{k). 


(26) 


(27) 


(28) 


Example 3. The long-range interaction 

K2{n,m) = - - 

\n — m\°‘ 

with the non-integer parameter a > 3, gives (see Theorem 8.7 in [5]) the relation 

K2(k) - K2(0) = -C(a - 2) fc" + ... , 


(29) 


( 30 ) 
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where is the Riemann zeta-function. 


Proposition 2 allows to demonstrate the close relation between the lattice structure and the 
fractional gradient non-local continuum. Let us describe the well-known special cases. 

Lattice equations fl2^ have two parameters g 2 and ga- The corresponding continuum equa¬ 
tion fl25D have two finite parameters C 2 and If we use 

5'2 = 4iL, 5'a = 0. 


then 

and equation 

then 


C 2 = Cl = Kd^M, = 0, 
gives equation flTT]) . If we assume that 

5-2 = 5'a = 4 iP, 


and we get relation 

d‘^u{x, t) 


Co = Cf = 


= cl Au{x, t) 


Kd^ 


On' 


2 ja-2 


r(a 1 ) 


9 ga-2 fi2 1 


r(cr -|- 1) 


p 


(31) 


(32) 


where Cg = \JE/p is the elastic bar velocity. Let us give a remark about the scale parameter 
Is{ol). Equation fIMl) can lead to incorrect conclusion about the behavior of the scale parameter 


Izia) = 


Ca P Co 


E Cl 

for d —0 in the case 0 < a < 2. Using Cl = Kd‘^/M, the parameter 

2Kd^ 


Kic^) = 


T{a + 1)CIM' 


(33) 

can be written as 

(34) 


Using that the value of Cl is finite, then behavior of the parameter ll{a) for d —)■ 0 has the 
same form for a > 2 and 0 < a < 2, such that ll{a) is proportion to d“. Therefore we assume 
that the range of validity of alpha parameter is arbitrary real positive a. 

For a = 4 equation fl3^ is the usual equations of the gradient elasticity models 


( 35 ) 


d‘^u{x,t) ^2 A / 3 a 2 / 3 

—— = Cg Au{x, t) - — A u{x, t) + -/(x). 

The correspondent stress-strain relation for linear one-dimensional elasticity has the form 

a{x,t) = e(^ 1 — llA^e{x,t), 











where a{x,t) is the stress, e{x,t) is the strain, and Is is the scale parameter. 

In general, the coupling constants g 2 and Qa are independent. Therefore the sign of the 
coupling constant Qq (including the case a = 4) may differ from the sign of the constant 
g 2 = 4: K. For a = 4 the second-gradient parameter is defined by the relation 

2 ^ \g4\ (P 
* 48 iF ’ 

where the sign in front of the factor is determined by the sign of the coupling constant g^. If 
the constant g^ is positive then we get the gradient elasticity model with negative sign |39] . 

As a result the second-gradient model with negative and positive sign can be derived from 
a microstructure of lattice particles by suggested approach. The suggested approach as shown 
above uniquely leads to second-order strain gradient terms that are preceded by the positive 
and negative signs. 



5 Stationary Solution for Fractional Gradient Elasticity 

We can consider more general model of lattice with long-range interaction, where all particles 
are displaced from its equilibrium in one direction, and the displacement of particles is described 
by a scalar held M(r,f), where r G M"' (n = 1,2,3). The correspondent continuum equation of 
the fractional elasticity model is 

^ ^ 2 ’ = C 2 AM(r, t) - Ca (-A)“/2y(r, t) + -/(r), (37) 

at^ p 

where r and r = |r| are dimensionless variables. 

Let us consider the static case {du{r, t)/dt = 0, i.e. M(r, t) = M(r)) in this fractional elasticity 
model. Then equation (l37ll has the form 

C 2 AM(r) - Ca (-A)"/^M(r) -h -/(r) = 0. (38) 

P 

We can use the Fourier method to solve fractional differential equation fl38|) . which is based 
on the relations 

J^[(-A)“/2^,(r)](k) = |k|“h(k), J’[An(r)](k) =-k2h(k). (39) 

Applying the Fourier transform PF to both sides of (l38il and using ([39]), we have 

(.FM)(k) = - [C 2 |k|2 + |k|“^)-' (-F/)(k). (40) 

P 

Equation fl38D (see, for example. Section 5.5.1. in [H]) has a particular solution that can 
be represented in the form of the convolution of the functions G”(|r|) and /(|r|) as follow 

«(r) = i [ G';(r-r')/(r')<iV, (41) 

P Jmx 
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where G^{r) is the Green function (see Section 5.5.1. in [H]) of the form 


G'lir) = T- 


(Gs |k |2 + |kr) 


-1 


(G2|kp + G„|k|“) ^ e+'(^’")ci’"k. (42) 


We can use the relation 

r Ir) poo 

/ e-<‘‘''>/(|k|)<i”k = )|G^ / /(A)A"/V„/ 2 -i(A|r|)<iA (43) 

jR’^ Jo 

that holds (see Lemma 25.1 of [13]) for any suitable function / such that the integral in the 
right-hand side of fl43l) is convergent. Here is the Bessel function of the hrst kind. 

Using relation fl43|l . the Green function fl4^ can be represented (see Theorem 5.22 in |14] i 
in the form of the integral with respect to one parameter A as 


r|( 2 -n )/2 ro- Xn/2 d\ 

(27r)G2 I G2A2 + G„A“ ’ 


(44) 


where n = 1, 2, 3, and J(n- 2)/2 is the Bessel function of the hrst kind. 
For the 3-dimensional case (n = 3), we can use 



Ji/ 2 iz) = \ — sin(^), 

" TIZ 


and we have 


G^fr) = 


A sin(A|r|) dX 


27r2|r|7o G2A2 + G„A^ 

For the 1-dimensional case (n = 1), we can use 



J-l/2{z) = \ — cos(z), 


TTZ 


and we have (see Theorem 5.24 in [Tl|) the function 

cos(A|r|) dX 


1 

Gl(r) = - 


7rJ„ C,A2 + C„A‘ 


(45) 

(46) 


(47) 


(48) 


Let us determine the deformation of an inhnite elastic continuum, when a force is applied 
to a small region of the medium. This is the well-known Thomson’s problem [10|. We solve 
Thomson’s problem in the framework of the fractional elasticity model. If we consider the 
deformation for |r|, which are larger compare with the size of the region, we can suppose that 
the force is applied at a point. In this case, we have 


/(r) = /o5(r) = fo6{x)6{y)6{z). 


(49) 
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Then the displacement held u{r) of fractional elasticity has a simple form of the particular 
solution that is proportional to the Green function 


M(r) = ^Gl{r). 

Therefore, the displacement held for the case 0491) has the form 

. . r Asin(A|r|) 

2n^p\v\]^ + 


(50) 


u\r} ^ 


d\. 


(51) 


The asymptotic behavior |r| —)■ cx) of the displacement held M(r) in the model described by 
(El]) with fl49|) . is given by 


u{r) ^ 


/o r(2 — a) sin(7rQ;/2) 


1 

J.|3 —Q 


(a<2), 


(52) 


M(r) ^ 


P |r 


(a > 2). 


(53) 


Note that the asymptotic behavior |r| —>■ cx) does not depend on the parameter a for a > 2. 
In the case a <2 the displacement held on the long distances is determined only by term with 
the fractional Laplacian of the order a. 

The asymptotic behavior |r| —)■ 0 of the displacement held M(r) that is described by equation 
cm. where the force /(r) is applied at a point fH9|) . is given by 


M(r) ^ 
u{r) 


r \ 1 /o 

“W ~ 0 2 11 

27r"' p |r| 

/or((3-«)/2) 

' pCaT{a/2) 

fo 


27iapc^p sin(37r/Q;) 


(a < 2), 

(54) 

3-a’ (2 < CK < 3), 

(55) 

, , , , (a > 3). 

(56) 

(Sn/a) 


Here the Euler’s rehection formula for Gamma function is used. Note that the asymptotic 
behavior |r| —)■ 0 does not depend on the parameter a for a < 2. In the case a > 2, the dis¬ 
placement held on the short distances is determined only by term with the fractional Laplacian 
of the order a. 

The functions 


u{x) = 


A sin(Aa:) 


dX 


for the diherent orders of 1 < a 


(T Jo C 2 A2 + A- 
< 6 and with C 2 = Ca = ^ are present on Figures 1-4, 


where x = |r|. Figures 1,2,4 allows us to see that the held u{x) tends to a constant value 
(m(x) —)■ const) at a; ^ 0 for the parameters a > 3 (a = 3.6, a = 4.1, a = 5.2 and a = 5.6). 
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Figures 2 and 3 demonstrate that the asymptotic behavior of the type u{x) ~ for the 

held u{x) at X —)■ 0 for the parameters 2 < a < 3 (a = 2.6, a = 2.7). Figures 3 and 4 show that 
the asymptotic behavior of the type uix) ~ 1/x for the held uix) at x —)■ 0 for the parameters 
0 < a < 2 (a = 1.4, a = 1.9). 

We can determine the deformation of an inhnite non-local elastic continuum, when a pair 
of forces with equal in magnitude and oppositely directed is applied to a small region of the 
medium. We assume that these forces are separated by small distance. If we consider the 
deformation for |r|, which are larger compare with the size of the region, we can suppose that 
two force are applied at two points such that 

/(r) =/o^(r + a)-/o5(r-a). (57) 


This problem is analogous to a dipole system as a pair of electric charges of equal magnitude 
but opposite sign, separated by small distance. Then the displacement held u{r) of fractional 
elasticity has a form of the particular solution 

w(r) = (^(((r + a)-^"(r-a)). (58) 

Therefore, the displacement held for the case fl57D has the form 


M(r) 


1 /o 

POO 

1 

^sin(A r -|- a ) sin(A r — a ) 

27r2 p J 

1 

r -I- a 


r — a 


A 

^2 A2 + A“ 


dX. 


(59) 


For |r| S> |a|, we can use |r — a| — |r -|- a| 2 |a| cos6', where 6 is the angle between r and a, 
and |r — a| |r -|- a| |rp. Then the displacement held can be represented in the form 


/o|a|cos6 ' A sin(A|r|) /o A sin(A |a| cos0) cos(A|r|) 

vrVIrp Jo C' 2 A 2 -kCoA" tt^ p\r\ Jq (^2 A^Co A“ 


For the suggested lattice equation and correspondent continuum limit, we can use all pos¬ 
itive values of alpha parameter. There is no reason to limit of the range alpha values in the 
used fractional diherential equations for the suggested form of fractional gradient and integral 
elasticity. It allows us to state that the range of validity of alpha parameter is arbitrary real 
positive values. As a result, we can distinguish two following particular cases in the fractional 
elasticity model described by fl38|) : (1) fractional integral elasticity (a < 2); (2) fractional gra¬ 
dient elasticity {a > 2). Note that for the hrst case the order of the fractional Laplacian in 
equation fl37)l is less than the order of the term related to the usual Hooke’s law. In the second 
case the order of the fractional Laplacian is greater of the order of the term related to the 
Hooke’s law. 
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6 Conclusion 


We suggest lattice models with long-range interaction of power-law type as microscopic model of 
fractional non-local elastic continuum. The continuum equations with derivatives of non-integer 
orders are directly derived from the suggested lattice models. We prove that the fractional 
gradient and fractional integral models can be derived from lattice models with long-range 
particle interactions. The suggested approach uniquely leads to second-order and fractional- 
order strain gradient terms that are preceded by the positive and negative signs. Fractional 
calculus allows us to obtain exact analytical solutions of the fractional differential equations for 
continuum models of a wide class of material with fractional gradient and fractional integral 
non-locality. A characteristic feature of the behavior of a fractional non-local continuum is the 
spatial power-tails of non-integer orders in the asymptotic behavior. The fractional elasticity 
models, which are suggested in this paper to describe complex materials with fractional non- 
locality, can be characterized by a common or universal spatial behavior of elastic materials 
by analogy with the universal temporal behavior of low-loss dielectrics im HU uni HZ]. The 
asymptotic behavior fl52p and 0541) allows us to state that fractional integral elasticity effects 
are important on the macroscopic scales. The asymptotic behavior (15^ and 05511561) allows us 
to state that fractional gradient elasticity effects are very important for the mesoscopic and 
nano scales. As a results the fractional gradient elasticity models can be very important for 
nanomechanics [HI 091001 011 02] of nonlocal materials with long-range particle interactions. 
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Appendix 1: Proof of Proposition 1 

To derive the equation for the field u{k,f), we multiply equation ffTOj) by exp(—zfcnd), and 
summing over n from — oo to J-oo. Then 

H-oo 72 -\-oo +00 

n=—oo n=—oo n=—oo 


The first term on the right-hand side of 0611) is 

H-OO H-OO 

K- Y e-^'^^’^K2{n,m)un = K ■ Y 


^n-|-l 2'U7 t, -|- lifi—i I 


H-oo 


X 


+00 


H-OO 




+ K. Y, 
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H-oo +00 +00 

m=—oo n=—oo j=—oo 

Using the definition of u{k,t), we obtain 
+00 

K ■ J2 m)un = K ■ t) - 2u{k, t) + e-*“M(A;, t)) = 

n=—oo 

= K ■ + e-*“ - 2^u{k, t) = 2K ■ (^cos {kd) - l) u{k, t). 

As a result, we have 

+c»o / JcrJ\ 

K- e~''^‘^'^K 2 {n^m)un =—^K ■ sw? i — \ u{k^t). (62) 

n=—oo ^ ^ 

Substitution of (16^ into (ET]) gives 

(^y) (“n(t))}- (63) 

For d —)■ 0, the asymptotic behavior of the sine is sin{kd/2) = kd/2 + 0{{kd)^). Then 

—4 sin^ ~ {kd)"^ + 0{{kdY). 

Using the finite parameter = Kd?/M^ the transition to the limit d —)• 0 in equation fl6^ 
gives 

^ = -ClPu{k, t) + ^J^{F{x)}, (64) 

where is defined by fll2l) . The inverse Fourier transform of fl6T)) has the form 

^ _c^2^-i{fc2a(<:,i)}+ )/(x), (65) 

where f{x) = F{x)/{A d) is the force density, and p = M/{Ad) is the mass density. Then using 

F~^{u{k,t)} = u{x,t), F~^{k'^u{k,t)} = —Au{x,t), 

we obtain the continuum equation ffTTj) . This ends the proof. 
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Appendix 2: Proof of Proposition 2 

To derive the equation for the held u{k,t), we multiply equation fl2^ by exp{—iknd), and 
summing over n from — oo to +oo. Then 


H-CXD 




+CXD +00 




n=—oo m=—oo 
m^n 


+ 00 +00 


9a ^ ^ e -m) Un - Mm j + 


+ 00 


—iknd 


F{n). 


n=—oo m=—oo 
m^n 


The left-hand side of fEB]) gives 


( 66 ) 


+ 00 

E 


.-ikndd'^Mt) _ g _ d'^u{k,t) 


dt2 Qt2 ^ Qt2 ’ 

n=—oo n=—oo 

where u{k,t) is dehned by (jS]). The second term of the right-hand side of fl66|) is 


(67) 


+00 


e-“”‘‘F(7i) = Ji{F(n)}. 


( 68 ) 


The limit for the hrst term on the right-hand side of fl66|) is described in Proposition 1. 
The second term on the right-hand side of (1661) with a multiplier is 


+00 +00 

E E 

n=—oo m=—oo 
m^n 


e - m) ( u„. - w™,) = 


+00 +00 


+ 00 +00 




n=—oo m=—oo 
m^n 


n=—oo m=—oo 
m^n 


Using (|3]), the hrst term in the right-hand side of (l69|) gives 


(69) 


e-“”V 


+ 00 +CXD +00 +00 

^ ^ e-^^'^'^K^{n-m)un= Y. 

n=—oo m=—oo n=—oo m'=—oo 

m^n mVO 

where we use (IT^ and Ka{m' + n — n) = Kaijn'), and 

+00 


= i^a(0)M(fc,t), (70) 


ko,{kd)= Y e-‘'="''ir„(n)= J-A{i^a(n)}. 


ri.=—OO 
n/0 


(71) 
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Note that 


+00 +00 H-OO +00 

^ ^ - m) = 


n= — 00 772= —oo 
m^n 


n=—oo m=—oo 

n^m 


+ 00 +00 

= e-“i^aK) u^e-^’^^^ = ka{kd)u{k,t), (72) 

n'=—oo m=—oo 

nVO 

where Ka{m — {n' + m)) = Kaiji') is used. 

As a result, equation (1661) has the form 

= 92 (i^2(0) - k^{kd)) u{Kt)+g. (i^a(O) - k^{kd)) u{k,t)+:F^{F{n)}, (73) 

where Fi^{F{n)} is an operator notation for the Fourier series transform of F{n). 

The Fourier series transform Fa of fl22|) gives fl73l) . We will consider the limit d —)■ 0. Using 
f[T6|l . equation fl73|l can be written as 

Fu(k, () = ^ %A(k) n{k, t) + ^ t„,i(fc) u(k, t} + T_f^{F(,!)}, (74) 

where we use (120|) . the Proposition 1 for K 2 {n,m), and the following notations 


Ta,A{k) = 


\k\^ + d^ 0{\k\^+^), 


2F(a + l) 
f2Aik) = -k^ + d‘^0{k). 


In the limit d —)■ 0, we get 


fkk)=C%Aik) = 


2F(a + l) 
f2(k) = C%A(k) = 


\k\“, 


(75) 

(76) 

(77) 

(78) 


As a result, equation fITil) in the limit d —)■ 0 gives 


where u{k,t) 
of d79|) is 


r)^ n n 1 

—U(k, () = ^ %(k) u(k, t) + ^ t(k) nik, t) + (79) 

= Cu{k,t) and F{k) = F{F{x)} = CFA{F{n)}. The inverse Fourier transform 


r X 92 d^ _, , , , Oad" ^ . 1 ^ / X 

ukF) + k{x) u{x,t) + -/(^), 


(80) 
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where f{x) = F{x)/{Ad) is the force density, the operators 72{x) and Ta{x) are dehned by 


%{x) = = A, r.(i) = T-^{t{k)} = (81) 

Here, we nse the connection between the Riesz’s fractional Laplacian and its Fonrier transform 
[13] in the form 

F[{-A)^/\{x)]{k) = \k\^u{k). (82) 

Using the finite parameters C 2 and Cq,, which are defined by (1211) . the snbstitntion of flST]) into 
fl5U]l gives continunm eqnations (1211) . This ends the proof. 
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